Quantum Spinodal Decomposition in Multicomponent Bose-Einstein Condensates 
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We investigate analytically the non-equilibrium spatial phase segregation process of a mixture of 
alkali Bose-Einstein condensates. Two stages (I and II) are found in analogy to the classical spinodal 
decomposition. The coupled non-linear Schrodinger equations enable us to give a quantitative picture 
of the present dynamical process in a square well trap. Both time and length scales in the stage I 
are obtained. We further propose that the stage II is dominated by the Josephson effect between 
different domains of same condensate different from scenarios in the classical spinodal decomposition. 
Its time scale is estimated. 
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The recent realizations of two |^J^ and three y] 
component alkali Bose-Einstein condensates (EEC's) in 
one trap provide us with new systems to explore the 
physics in otherwise unachievable parameter regimes 
[H-p|. These systems have two noticeable advantages: 
the easy control of experimental parameters and the rel- 
ative simplicity of the mathematical description. A direct 
comparison between theoretical calculations and exper- 
imental observations can be made. The purpose of the 
present paper is to explore one of main questions in the 
non-equilibrium statistical physics using those new sys- 
tems: spinodal decomposition in a binary-solution sys- 
tem. It is a typical example of phase ordering dynamics: 
the growth of order through domain coarsening when a 
system is quenched from the homogeneous phase into a 
broken-symmetry phase H. Systems quenched from a 
disordered phase into an ordered phase do not order in- 
stantaneously. Instead, different length scales set in as 
the domains form and grows with time, and different bro- 
ken symmetry phases compete to select the equilibrium 
state. We show that it is possible to have an analogous 
spinodal decomposition in EEC's, which manifests all 
the main phenomenology except those constrained by the 
trap size. In the following, we shall study the main fea- 
tures in this dynamical evolution process, starting from 
the homogeneous unstable state of EEC's. To differen- 
tiate the present situation from the usual one, we shall 
call the present one the quantum spinodal decomposition, 
and the previous ones the classical spinodal decomposi- 
tion. 

We start from the time dependent non-linear 
Schrodinger equations 
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Here i^j{x, i), toj, Uj with j — 1, 2 are the effective wave 
function, the mass, and the trapping potential of the 
jth condensate. The repulsive interaction between the 
jth condensate atoms is specified by Gjj, and that be- 
tween 1 and 2 by G12. The Lagrangian multipliers, the 
chemical potentials /ii and /i2, are fixed by the relations 
j d?'x\%ljj{x,t)\'^ = Nj ,j — 1,2, with Nj the number of 
the jth condensate atoms. Eq. (1) and (2) are mean field 
equations, since we treat the effective wave functions as c 
numbers, corresponding to the Hartree-Fock-Bogoliubov 
approximation. They provide a good description for the 
slow dynamics in both alkali EEC's and the superfluid 
Hc4 pi on length scale larger than the range of micro- 
scopic interaction. 

Experimentally, the trapping potentials {Uj} are sim- 
ple harmonic in nature. For the sake of simplicity and 
to illustrate the physics we shall consider a square well 
trapping potential Uj = U: zero inside and large (infi- 
nite) outside, unless otherwise explicitly specified. 

We consider the strong mutual repulsive interaction 
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In this regime the equilibrium state for two Bose-Einstein 
condensates is a spatial segregation of two condensates, 
where two phases, the weakly and strong segregated 
phases, characterized by the healing length and the pen- 
etration depth, have been predicted 0. We shall use Eq. 



(1) and (2) under the condition (3) to study a highly 
non-hnear dynamical process: The two condensates are 
initially in a homogeneously mixed state, then eventu- 
ally approach the phase segregated state. In the same 
mean field treatment as in Ref. [B , we find that this dy- 
namical process can be classified into two main stages: 
The initial highly non-equilibrium dynamical growth in 
the stage I, where the dynamics is governed by the fastest 
growth mode, and the stage II of approach to equilibrium 
where the dynamics is governed by the slowest mode. 
The stage II is typical of a relaxation process near equi- 
librium. However, we shall show again it is governed by 
a quantum effect, namely, the Josephson effect. 

Stage I: Fastest Growth Mode. The coupled non- 
linear Schrodinger equations have an obvious homoge- 
neous solution: Inside the trap the condensate densities 
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Pjo-i Pjo = Nj/V , with V the volume of the 



square well potential trap, and the chemical potentials 

Ml = GiiPw + G12P20 and fi2 = G22P20 + Gi2Pio- This 
is the initial condition of the present problem. 

To look for the fastest growth mode out of the ho- 
mogeneous state, we start with small fluctuations from 
the homogeneous state. This is consistent with the usual 
stability analysis Q. Our approach here is to emphasize 
the connection with the physics of the classical spinodal 
decomposition and the role played by the Josephson re- 
lationships. Define 
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and define the density fluctuations Spj ^ pj — pjo and 
the phase fluctuations 9j, and assume they are small: 
\Spj\/pj , |6'j| << 1. The definition of the phase fluctua- 
tions here has made use of the implicit assumption that 
there is no net current in the condensate. To the lin- 
ear order, we have from Eqs.(l,2), after eliminating the 
phase variables. 
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We look for the solution 
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with A, B constants. There are two branches of solu- 
tion for Eq. (5). For one branch, the frequency is al- 
ways real. For another branch which we denote by Ci;„, 
it can be imaginary. An imaginary frequency a;_ shows 
that the initial homogeneously mixed state is unstable. 



It is straightforward to verify that the sufficient condi- 
tion for the appearance of imaginary frequency for small 
enough wave number q is the validity of the inequal- 
ity Eq.(3). The modes with imaginary frequencies will 
then grow exponentially with time. Unlike the usual sit- 
uation near equilibrium, this growth from the present 
non-equilibrium homogeneously mixed state will be dom- 
inated by the fastest growth mode. This is precisely the 
same case as in the initial stage of the classical spinodal 
decomposition. 

To get a concrete understanding of the physical im- 
plications of Eq. (5) in the stage I we consider a case 
relevant to recent experiments where particles of the two 
condensates have the same mass, toi = TO2 = rn. In this 
case we find the wavenumber corresponding to the most 



negative uj 


2 is 
— max -^"^ 








2 "^ 

^max ~r2 


Jbl + 4pio/?2oG'llG22 


f GI2 
\G11G22 


-0 


-60 

(6) 



and 



■ ^ 2 



(7) 



with 60 = PioGii + p2oG22- The physics implied in Eqs. 
(6,7) is as follows. Starting from the initial homogeneous 
mixture of the two condensates, on the time scale given 

by 



tl = l/\tO- 
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domain patterns of the phase segregation with the char- 
acteristic length 



1 1 = l/<7r, 
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will appear. Particularly, for the weakly segregated 
phase of GI2/G11G22 — 1 — > 0, we have the length 
scale IJ^ — \f2l \J !s\ -t- A2, and for the strongly segre- 
gated phase of G\2l G\\G22 — 1 -^ 00, we have l~[^ = 
V2/VA TX;. Here A, ^ ^j/VGi2/VGiiG22 ~ 1 and 
^j — yTi /2m j I/pjqGjj are the penetration and heal- 
ing lengths in the binary EEC mixture Q . Those length 
and times scales can be measured experimentally. We 
will come back to the experimental situation below. 

After the stage I of fast growth into the domain pat- 
tern characterized by the length scale Ij, the system will 
gradually approach the true ground state of the complete 
phase segregation: one condensate in one region and the 
second condensate in another region. This stage is slow. 



dominated by the slowest mode, and is the subject of 
stage II. 

It is now evident that the stage I of the growth of bi- 
nary BEC's shares the same phenomenology of the initial 
stage of the classical spinodal decomposition: the dom- 
ination of the fastest growth mode, the appearance of 
domains of segregated phases, and the conservation of 
particle numbers. There are, however, two important 
differences. First, the dynamical evolution of the binary 
BEC's is governed by a coupled nonlinear time depen- 
dent Schrodinger equations, not by a nonlinear diffusion 
equation supplemented with the continuity equation, the 
Cahn-Hilliard equation. There is no external relaxation 
process for the present wave functions. Secondly, the en- 
ergy of the binary BEC's is conserved during the growth 
process, not as in the case of classical spinodal decompo- 
sition where the system energy always decreases. 

Stage 11: Merging and Oscillating between Domains. 
The BEC binary mixture occurs in a trap. This finite 
size effect of the droplet leads to the broken symmetry of 
the condensate profiles M, which tends to separate the 
condensates in mutually isolated regions. This implies 
that there is no contact between different domains of the 
same condensates formed in the stage I. The classical 
spinodal process involves diffusion. An estimate of the 
diffusion constant for the BEC system can be made from 
kinetics theory. The ratio of the time scales for quantum 
and classical particle transport is of the order of the ratio 
of the BEC cloud size to the de Broglie wavelength. This 
is much smaller than one for the experimental systems 
of interest. The classical diffusion process is thus not 
important. Because the domains are not connected and 
because the diffusion for the BEC mixture is extremely 
low, all the mechanisms for the late stage classical spin- 
odal decomposition are not applicable. We propose that 
it is the Josephson effect that is responsible for the ap- 
proach to equilibrium in the stage II. Two models for 
the Josephson effect, the 'rigid pendulum' model and the 
'soft pendulum' model will be discussed. They both give 
the same time scale when the 'Rabi' frequency is small. 

Let us consider the specific case of two domains of con- 
densate 1 separated by a domain with width d of con- 
densate 2. The ability of condensate 1 to tunnel through 
condensate 2 is described by the penetration depth A, as 
discussed in Ref. S] . Hence the probability of condensate 
1 to tunnel through condensate 2 can be estimated as 
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when p is much smaller than 1. The finiteness, though 
small, of the tunneling probability suggests that it is the 
Josephson effect responsible for the relaxation process in 



the stage II. The Josephson effect leads to the merging 
of two domains of the same condensate. The dynamics 
of such motion may be governed by the 'rigid pendulum' 
Hamiltonian for a Josephson junction pO[|: 
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where Ej is the Josephson coupling energy determined 
by the tunneling probability, n = (jix — ny)/2 is the par- 
ticle number difference between the numbers of particles, 
n^ and Uy, in the two domains, and Ec = dfi/dn is the 
'capacitive' energy due to interactions. In the absence of 
external constraints, fi = Equ. The phase difference 
between the two domains is conjugated to n, as in usual 
Josephson junctions. Under the appropriate condition, 
such as low temperature and smallncss of the capacitive 
energy, there may be an oscillation between the two do- 
mains of condensate 1 separated by condensate 2. In 
such a case, we may estimate that the oscillation period 
tii = 2t: Jujjp, with the so-called Josephson plasma fre- 
quency 1 10 
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For small tunneling probability, the Josephson junc- 
tion energy may be estimated as [lO[ Ej = 



1/3 
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g {nr/i) ^'^ (15aii/ao)'^'* ?iwo- Here wq is the harmonic 
oscillator frequency for condensate 1 in a harmonic trap, 
Co — y^h/miujQ is the corresponding oscillator length, 
the flu is the scattering length of condensate 1, and 
nx — n^ + Uy is the total number of particles in do- 
main X and y. The oscillatory time scale between the 
domains determines by the Josephson plasma frequency 



// 



2ao 



15aii nr 



2/15 



Wo 

2tt 



-d/A 



(12) 



The rigid pendulum Hamiltonian would give a good de- 
scription when n << ut- Another description of the 
Josephson effect uses the 'soft pendulum' Hamiltonian 
proposed in Ref. [|ll| . In the parameter regime relevant 
current experimental situations both of them give the 
same order of magnitude estimate for t//. 

Given the Josephson effect is the dominant mechanism 
in the stage II, the time scale to arrive at the ground state 
will be determined by the Josephson effect at the final 
two domains, in which the two domains of condensate 1 
is separated by the domain of condensate 2, in the man- 
ner of 1212 spatial configuration for the case of equal 
numbers of the condensates. The width d of each domain 



is then D/4, with D the size of the trap. According to 
the above analysis, the largest time scale determined by 
Eq.(12), the slowest mode, in the stage II is i//. The 
arrangement of 1 2 1 spatial configuration may also occur 
here, in which it is more likely for condensate 2 to tunnel 
through 1 to the edge of the trap, because of the larger 
tunneling probability. We next turn our attention to the 
experimental situation. 

The first question is whether the quantum spinodal 
decomposition discussed above can happen or not. In 
terms of the atomic scattering lengths of condensate 



atoms Qjj, the interactions are G,-, = 47r?i a 



for 8'^Rb is about 50A. The 



The typical value of ajj 

typical density realized for the binary BEC mixture is 

about pio ~ 10^^/cTO'^. Hence the healing length is 

e = ^{fi'/2m){l/GjjPjo) = VV(8™«-Pjo) - 3000A. 
For the different hyperfine states of ®'^Rb, it is known now 
|g,|| that Gi2/y/GiiG22 > 1- Hence the ground state of 
the phase segregated phase can be realized. Also, that 
both the healing length and the penetration depth are 
smaller than the trap size Q validates the application 
of results obtained in the square well trapping potential 
to the harmonic one. Therefore, the quantum spinodal 
decomposition can happen. 

Experimentally, starting from the initially homoge- 
neous state, after a short period of time a domain pat- 
tern does appear. Then a damped oscillation between 
the domain pattern has been observed. Eventually the 
binary BEC mixture sets to the segregated phase Q . If 
we take Gi2/\/GiiG22 = 1-04, the penetration depth is 
A = ^/^/(h^J7^^l^^^ - 1.5/im. The length scale 
li determined by Eq.(9) in the stage I is 1.5^m, which 
is the same order of magnitude in comparison with ex- 
perimental data H. This length is also comparable to 
the domain wall width seen experimentally. The corre- 
sponding time scale tj (Eq. (8)) is then 6 ms, again the 
same order of magnitude, though both estimated li and 
ti are somewhat smaller. For stage II, d/A is of the order 
of Ud, the number of domains formed in stage I. For the 
experiment of interest to use, Ud ~ 2. If we assume that 
the damped oscillation to equilibrium observed experi- 
mentally is the stage II discussed here, taking the total 
number of particle A^i = 10^, we find the period accord- 
ing to Eq. (12) is 30 ms, comparable to the experimental 
value. 

For the stage I, our above analysis shows that it is in- 
sensitive to the damping. At this moment we do not have 
a reliable estimation of the damping in the Josephson os- 
cillation, whose existence has been indicated experimen- 
tally. Nevertheless, given the uncertainty in the value of 



Gi2 or ai2, we conclude that the stage II of the quantum 
spinodal decomposition may have been observed. 

Periodic-like structures have also been observed in the 
phase segregation of spin 1 Na mixtures |^. For the 
definiteness of analysis we have presently focused on the 
binary case, but we believe a similar physical picture of 
quantum spinodal decomposition applies to that case as 
well. 

Finally, we point out that the length and time scales 
for Rb mixtures is controlled by the difference of r = 
G12/VG11G22 and 1. Since r is close to 1 experimentally, 
the data can provide for a very sensitive estimate of r- 
1. Since this parameter determines both the time and 
length scale of stage I, it is a self-consistent check on the 
physical picture provided here. 
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